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METHOD AND APPARATUS FOR FINE FREQUENCY SYNCHRONIZATION IN ~ 



MULTI-CARRIER DEMODULATION SYSTEMS 



FIELD OF THE INVENTION 



The present invention relates to methods and apparatus for 
10 performing a fine frequency synchronization in multi -carrier 
demodulation systems, and in particular to methods and appa- 
ratus for performing a fine frequency synchronization com- 
pensating for a carrier frequency deviation from an oscilla- 
tor frequency in a multi -carrier demodulation system of the 
15 type capable of carrying out a differential phase decoding 
of multi -carrier modulated signals, wherein the signals com- 
prise a plurality of symbols, each symbol being defined by 
phase differences between simultaneous carriers having dif- 
ferent frequencies. 

20 

BACKGROUND OF THE INVENTION 

In a multi carrier transmission system (MCM, OFDM), the ef- 

25 feet of a carrier frequency offset is substantially more 
considerable than in a single carrier transmission system. 
MCM is more sensitive to phase noise and frequency offset 
which occurs as amplitude distortion and inter carrier in- 
terference (ICI) . The inter carrier interference has the ef- 

30 feet that the subcarriers are no longer orthogonal in rela- 
tion to each other. Frequency offsets occur after power on 
or also later due to frequency deviation of the oscillators 
used for downconversion into baseband. Typical accuracies 
for the frequency of a free running oscillator are about ±50 

35 ppm of the carrier frequency, with a carrier frequency in 
the S-band of 2.34 Ghz, for example, there will be a maximum 
local oscillator (LO) frequency deviation of above 100 kHz 
(117.25 kHz). The above named effects result in high re- 
quirements on the algorithm used for frequency offset cor- 

40 recti on. 




DESCRIPTION OF PRIOR ART 

5 Most prior art algorithms for frequency synchronization di- 
vide frequency correction into two stages. In the first 
stage, a coarse synchronization is performed. In the second 
stage, a fine correction can be achieved. A frequently used 
algorithm for coarse synchronization of the carrier fre- 
10 quency uses a synchronization symbol which has a special 
spectral pattern in the frequency domain. Such a synchroni- 
zation symbol is, for example, a CAZAC sequence (CAZAC = 
^ Constant Amplitude Zero Autocorrelation). Through compari- 

2 son, i.e. the correlation, of the power spectrum of the re- 

CP 15 ceived signal with that of the transmitted signal, the fre- 
quency carrier offset can be coarsely estimated. These prior 
ru art algorithms all work in the frequency domain. Reference 

2? is made, for example, to Ferdinand Cl alien, Heinrich Meyr, 

"Synchronization Algorithms for an OFDM system for Mobile 
[7 20 Communication", ITG-Fachtagung 130, Codierung fur Quelle, 
m Kanal und ubertragung, pp. 105 - 113, Oct. 26-28, 1994; and 

Timothy M. Schmidl, Donald C. Cox, "Low-Overhead, Low- 
g Complexity [Burst] Synchronization for OFDM", in Proceedings 

of the IEEE International Conference on Communication ICC 
25 1996, pp. 1301-1306 (1996). 

For the coarse synchronization of the carrier frequency, 
Paul H. Moose, "A Technique for Orthogonal Frequency Divi- 
sion Multiplexing Frequency offset Correction", IEEE Trans- 

30 act i on On Communi cati ons , Vol . 42 , No . 10 , Octobe r 1994 , 
suggest increasing the spacing between the subcarriers such 
that the subcarrier distance is greater than the maximum 
frequency difference between the received and transmitted 
carriers. The subcarrier distance is increased by reducing 

35 the number of sample values which are transformed by the 
Fast Fourier Transform. This corresponds to a reduction of 
the number of sampling values which are transformed by the 
Fast Fourier Transform. 
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WO 9205646 A relates to methods for the reception of or- 
thogonal frequency division multiplexed signals comprising 
data which are preferably differentially coded in the direc- 
tion of the time axis. Phase drift of the demodulated sam- 
5 pies from one block to the next is used to indicate the de- 
gree of local oscillator frequency error. Phase drift is as- 
sessed by multiplying complex values by the complex conju- 
gate of an earlier sample demodulated from the same OFDM 
carrier and using the resulting measure to steer the local 
10 oscillator frequency via a frequency locked loop. 



_ SUMMARY OF THE INVENTION 

m 15 It is an object of the present invention to provide methods 
fl and apparatus for performing a fine frequency synchroniza- 

rU tion which allow a fine frequency synchronization compensat- 

2i ing for a carrier frequency deviation from an oscillator 

7 frequency in a MCM transmission system which makes use of 

H 20 MCM signals in which information is differential phase en- 
coded between simultaneous sub-carriers having different 
frequencies . 



in accordance with a first aspect, the present invention 
25 provides a method of performing a fine frequency synchroni- 
zation compensating for a carrier frequency deviation from 
an oscillator frequency in a multi-carrier demodulation sys- 
tem of the type capable of carrying out a differential phase 
decoding of multi -carrier modulated signals, the signals 
30 comprising a plurality of symbols, each symbol being defined 
by phase differences between simultaneous carriers having 
different frequencies, the method comprising the steps of: 

determining a phase difference between phases of the same 
35 carrier in different symbols; 



determining a frequency offset by eliminating phase shift 
uncertainties related to the transmitted information from 
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the phase difference making use of a M-PSK decision device; 
and 

performing a feedback correction of the carrier frequency 
5 deviation based on the determined frequency offset. 

In accordance with a second aspect, the present invention 
provides a method of performing a fine frequency synchroni- 
zation compensating for a carrier frequency deviation from 

10 an oscillator frequency in a multi -carrier demodulation sys- 
tem of the type capable of carrying out a differential phase 
decoding of multi -carrier modulated signals, the signals 
comprising a plurality of symbols, each symbol being defined 
by phase differences between simultaneous carriers having 

15 different frequencies, the method comprising the steps of: 

determining respective phases of the same carrier in differ- 
ent symbols; 

20 eliminating phase shift uncertainties related to the trans- 
mitted information from the phases to determine respective 
phase deviations making use of a M-PSK decision device; 

determining a frequency offset by determining a phase dif- 
25 ference between the phase deviations; and 

performing a feedback correction of said carrier frequency 
deviation based on the determined frequency offset. 

30 In accordance with a third aspect, the present invention 
provides an apparatus for performing a fine frequency syn- 
chronization compensating for a carrier frequency deviation 
from an oscillator frequency, for a multi -carrier demodula- 
tion system of the type capable of carrying out a differen- 

35 tial phase decoding of multi -carrier modulated signals, the 
signals comprising a plurality of symbols, each symbol being 
defined by phase differences between simultaneous carriers 
having different frequencies, the apparatus comprising: 
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means for determining a phase difference between phases of 
the same carrier in different symbols; 

M-PSK decision device for determining a frequency offset by 
5 eliminating phase shift uncertainties related to the trans- 
mitted information from the phase difference; and 

means for performing a feedback correction of the frequency 
deviation based on the determined frequency offset. 

10 

In accordance with a fourth aspect, the present invention 
provides an apparatus for performing a fine frequency syn- 
chronization compensating for a carrier frequency deviation 
from an oscillator frequency, for a multi-carrier demodula- 
15 tion system of the type capable of carrying out a differen- 
tial phase decoding of multi -carrier modulated signals, said 
signals comprising a plurality of symbols, each symbol being 
defined by phase differences between simultaneous carriers 
having different frequencies, the apparatus comprising: 

20 

means for determining respective phases of the same carrier 
in different symbols; 

M-PSK decision device for eliminating phase shift uncertain- 
25 ties related to the transmitted information from the phases 
to determine respective phase deviations; 

means for determining a frequency offset by determining a 
phase difference between the phase deviations; and 

30 

means for performing a feedback correction of the frequency 
deviation based on the determined frequency offset. 

The present invention relates to methods and apparatus for 
35 performing a fine frequency synchronization compensati ng for 
a carrier frequency deviation from an oscillator frequency. 
This fine frequency synchronization is preferably performed 
after completion of a coarse frequency synchronization, such 
that the frequency offsets after the coarse frequency syn- 
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chronization are smaller than half the sub-carrier distance 
in the MCM signal- Since the frequency offsets which are to 
be corrected by the inventive fine frequency synchronization 
methods and apparatus, a correction of the frequency offsets 
5 by using a phase rotation with differential decoding and de- 
mapping in the time axis can be used. The frequency offsets 
are detected by determining the frequency differences be- 
tween time contiguous sub-carrier symbols along the time 
axis. The frequency error is calculated by measuring the ro- 
10 tation of the I-Q cartesian coordinates of each sub-carrier 
and, in preferred embodiments, averaging them over all n 
sub-carriers of a MCM symbol. 

^ Firstly, the phase ambiguity or uncertainty is eliminated by 
m 15 using a M-PSK decision device and correlating the output of 
s 7j the decision device with the input signal for a respective 
ry sub-carrier symbol. Thus, the phase offset for a sub-carrier 
2* symbol is determined and can be used for restructuring the 
/ frequency error in form of a feed-backward structure. Alter- 
20 natively, the phase offsets of the sub-carrier symbols of 
one MCM symbol can be averaged over all of the active carri- 
es ers of a MCM symbol, wherein the averaged phase offset is 
^ used to restructure the frequency error. 

25 In accordance with the present invention, the determination 
of the frequency offset is performed in the frequency do- 
main. The feedback correction in accordance with the inven- 
tive fine frequency synchronization is performed in the time 
domain. To this end, a differential decoder in the time do- 

30 main is provided in order to detect frequency offsets of 
sub-carriers on the basis of the phases of timely successive 
sub-carrier symbols of different MCM symbols. 



35 BRIEF DESCRIPTION OF THE DRAWINGS 



In the following, preferred embodiments of the present 
vention will be explained in detail on the basis of 
drawings enclosed, in which: 



i n- 
the 



10 
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Figure 1 shows a schematic overview of a MCM transmis- 
sion system to which the present application 
can be applied; 

Figures 2A and 2b show schematic views representing a 
scheme for differential mapping in the time 
axis and a scheme for differential mapping in 
the frequency axis; 

Figure 3 shows a functional block diagram for performing 
a differential mapping in the frequency axis; 



Figure 4 shows a representation of time variation of all 
m 15 sub-carriers in MCM symbols; 

Li! 

ry Figure 5 shows a QPSK-constellation for each sub-carrier 

^ with a frequency offset; 

^ 20 Figure 6 shows a general block diagram illustrating the 

£1 position of the inventive fine frequency syn- 

y3 chronization device in a MCM receiver; 

ci 

Figure 7 shows a block diagram of the fine frequency er- 
25 ror detector shown in Figure 6; 

Figure 8 shows a block diagram of a MCM receiver com- 
prising a coarse frequency synchronization unit 
and a fine frequency synchronization unit; 

30 

Figure 9 shows a block diagram of a unit for performing 
a coarse frequency synchronization; 

Figure 10 shows a schematic view of a reference symbol 
35 used for performing a coarse frequency synchro- 

nization; 



Figure 11 shows a schematic view of a typical MCM signal 
having a frame structure; 
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shows scatter diagrams of the output of an dif- 
ferential d(i-mapper of a mcm receiver for il- 
lustrating i:he effect of an echo phase offset 
correction; ^ 

shows a schematic block diagram for illustrat- 
ing the position and the functionality of an 
echo phase offset correction unit; 

shows a schematic block diagram of a preferred 
form of an echo phase offset correction device; 
and 

shows schematic views for illustrating a pro- 
jection performed by another echo phase offset 
correction algorithm. 

20 DETAILED DESCRIPTION OF THE EMBODIMENTS 

Before discussing the present invention in detail, the mode 
of operation of a MCM transmission system is described re- 
ferring to figure 1. 

25 

Referring to Figure 1, at 100 a MCM transmitter is shown 
that substantially corresponds to a prior art MCM transmit- 
ter. A description of such a MCM transmitter can be found, 
for example, in William Y. zou, Yiyan wu, "COFDM: AN 
30 OVERVIEW", IEEE Transactions on Broadcasting, vol. 41, No. 
1, March 1995. 

A data source 102 provides a serial bitstream 104 to the MCM 
transmitter. The incoming serial bitstream 104 is applied to 
35 a bit-carrier mapper 106 which produces a sequence of spec- 
tra 108 from the incoming serial bitstream 104. An inverse 
fast Fourier transform (IFFT) 110 is performed on the se- 
quence of spectra 108 in order to produce a MCM time domain 
signal 112. The MCM time domain signal forms the useful MCM 



Figure 12 



10 



Figure 13 



Figure 14 



15 Figure 15 



10 



m 
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symbol of the MCM time signal. To avoid intersymbol inter- 
ference (ISI) caused by multipath distortion, a unit 114 is 
provided for inserting a guard interval of fixed length be- 
tween adjacent MCM symbols in time, in accordance with a 
preferred embodiment of the present invention, the last part 
of the useful MCM symbol is used as the guard interval by 
placing same in front of the useful symbol. The resulting 
MCM symbol is shown at 115 in Figure 1 and corresponds to a 
MCM symbol 160 depicted in Figure 11. 



Figure 11 shows the construction of a typical MCM signal 
having a frame structure. One frame of the MCM time signal 
is composed of a plurality of MCM symbols 160. Each MCM sym- 
bol 160 is formed by an useful symbol 162 and a guard inter- 
15 val 164 associated therewith. As shown in Figure 11, each 
frame comprises one reference symbol 166. The present inven- 
tion can advantageously be used with such a MCM signal, how- 
^ ever, such a signal structure being not necessary for per- 

r forming the present invention as long as the transmitted 

^ 20 signal comprises a useful portion and at least one reference 
hi symbol . 



in order to obtain the final frame structure shown in Figure 
11, a unit 116 for adding a reference symbol for each prede- 
2 5 termined number of MCM symbols is provided. 

In accordance with the present invention, the reference sym- 
bol is an amplitude modulated bit sequence. Thus, an ampli- 
tude modulation of a bit sequence is performed such that the 

30 envelope of the amplitude modulated bit sequence defines a 
reference pattern of the reference symbol. This reference 
pattern defined by the envelope of the amplitude modulated 
bit sequence has to be detected when receiving the MCM sig- 
nal at a MCM receiver, in a preferred embodiment of the pre- 

35 sent invention, a pseudo random bit sequence having good 
autocorrelation properties is used as the bit sequence that 
is amplitude modulated. 
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The choice of length and repetition rate of the reference 
symbol depends on the properties of the channel through 
which the MCM signal is transmitted, e.g. the coherence time 
of the channel- in addition, the repetition rate and the 
5 length of the reference symbol, in other words the number of 
useful symbols in each frame, depends on the receiver re- 
quirements concerning mean time for initial synchronization 
and mean time for resynchronization after synchronization 
loss due to a channel fade. 

10 

The resulting MCM signal having the structure shown at 118 
in Figure 1 is applied to the transmitter front end 120. 
Roughly speaking, at the transmitter front end 120, a digi- 
tal/analog conversion and an up-converting of the MCM signal 
15 is performed- Thereafter, the MCM signal is transmitted 
through a channel 122. 

Following, the mode of operation of a MCM receiver 130 is 
shortly described referring to Figure 1. The MCM signal is 
20 received at the receiver front end 132- In the receiver 
front end 132, the MCM signal is down-converted and, fur- 
thermore, an analog/digital conversion of the down-converted 
signal is performed. 

25 The down-converted MCM signal is provided to a symbol 
frame/carrier frequency synchronization unit 134. 

A first object of the symbol frame/carrier frequency syn- 
chronization unit 134 is to perform a frame synchronization 
30 on the basis of the amplitude-modulated reference symbol. 
This frame synchronization is performed on the basis of a 
correlation between the amplitude-demodulated reference sym- 
bol and a predetermined reference pattern stored in the MCM 
receiver . 

35 

A second object of the symbol frame/carrier frequency syn- 
chronization unit is to perform a coarse frequency synchro- 
nization of the MCM signal. To this end, the symbol 
frame/carrier frequency synchronization unit 134 serves as a 
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coarse frequency synchronization unit for determining a 
coarse frequency offset of the carrier frequence caused, for 
example, by a difference of the frequencies between the lo- 
cal oscillator of the transmitter and the local oscillator 
5 of the receiver. The determined frequency is used in order 
to perform a coarse frequency correction. The mode of opera- 
tion of the coarse frequency synchronization unit is de- 
scribed in detail referring to Figures 9 and 10 hereinafter. 

10 As described above, the frame synchronization unit 134 de- 
termines the location of the reference symbol in the MCM 
symbol. Based on the determination of the frame synchroniza- 
tion unit 134, a reference symbol extracting unit 136 ex- 
tracts the framing information, i.e. the reference symbol, 

15 from the MCM symbol coming from the receiver front end 132. 
After the extraction of the reference symbol, the MCM signal 
is applied to a guard interval removal unit 138. The result 
of the signal processing performed hereherto in the MCM re- 
ceiver are the useful mcm symbols. 

20 

The useful MCM symbols output from the guard interval re- 
moval unit 138 are provided to a fast Fourier transform unit 
140 in order to provide a sequence of spectra from the use- 
ful symbols. Thereafter, the sequence of spectra is provided 
25 to a carrier-bit mapper 142 in which the serial bitstream is 
recovered. This serial bitstream is provided to a data sink 
144. 

Next, referring to Figures 2A and 2B, two modes for differ- 
30 ential mapping are described, in Figure 2a, a first method 
of differential mapping along the time axis is shown. As can 
be seen from Figure 2A, a MCM symbol consists of K sub- 
carriers. The sub-carriers comprise different frequencies 
and are, in a preferred embodiment, equally spaced in the 
35 frequency axis direction. When using differential mapping 
along the time axis, one or more bits are encoded into phase 
and/or amplitude shifts between two sub-carriers of the same 
center frequency in adjacent MCM symbols. The arrows de- 
picted between the sub-carrier symbols correspond to infor- 
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mation encoded in amplitude and/or phase shifts between two 
sub-carrier symbols. 

A second method of differential mapping is shown in Figure 
2b. The present invention is adapted for MCM transmission 
system using the mapping scheme shown in Figure 2B. This 
mapping scheme is based Ion a differential mapping inside one 
MCM symbol along the frequency axis. A number of MCM symbols 
200 is shown in Figure KB. Each MCM symbol 200 comprises a 
10 number of sub-carrier symbols 202. The arrows 204 in Figure 
2b illustrate information encoded between two sub-carrier 
symbols 202. As can be seen from the arrows 204, this map- 
ping scheme is based onla differential mapping within one 
MCM symbol along the frequency axis direction, 
15 

In the embodiment shown in Figure 2B, the first sub-carrier 
(k=0) in an MCM symbol 200 is used as a reference sub- 
carrier 206 (shaded) such that information is encoded be- 
tween the reference sub-carrier and the first active carrier 
20 208. The other information of a MCM symbol 200 is encoded 
between active carriers, respectively. 



Thus, for every MCM symbol an absolute phase reference ex- 
ists. In accordance with Figure 2B, this absolute phase ref- 

25 erence is supplied by a reference symbol inserted into every 
MCM symbol (k=0) . The reference symbol can either have a 
constant phase for all MCM symbols or a phase that varies 
from MCM symbol to MCM symbol. A varying phase can be ob- 
tained by replicating the phase from the last subcarrier of 

30 the MCM symbol preceding in time. 

in Figure 3 a preferred embodiment of a device for perform- 
ing a differential mapping along the frequency axis is 
shown. Referring to Figure 3, assembly of MCM symbols in the 
35 frequency domain using differential mapping along the fre- 
quency axis according to the present invention is described. 

Figure 3 shows the assembly of one MCM symbol with the fol- 
lowing parameters: 
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NFFT designates the number of complex coefficients of the 
discrete Fourier transform, number of subcarriers re- 
spectively. 

5 

K designates the number of active carriers. The refer- 
ence carrier is not included in the count for K. 

According to Figure 3, a quadrature phase shift keying 
10 (QPSK) is used for mapping the bitstream onto the complex 
symbols. However, other M-ary mapping schemes (MPSK) like 2- 
PSK, 8-PSK, 16-QAM, 16-APSK, 64-APSK etc. are possible. 

Furthermore, for ease of filtering and minimization of 
15 aliasing effects some subcarriers are not used for encoding 
information in the device shown in Figure 3. These subcarri- 
ers, which are set to zero, constitute the so-called guard 
bands on the upper and lower edges of the MCM signal spec- 
trum. 

20 

At the input of the mapping device shown in Figure 3, com- 
plex signal pairs b0[k], bl[k] of an input bitstream are re- 
ceived. K complex signal pairs are assembled in order to 
form one MCM symbol. The signal pairs are encoded into the K 
25 differential phase shifts phi [k] needed for assembly of one 
MCM symbol. In this embodiment, mapping from Bits to the 0, 
90, 180 and 270 degrees phase shifts is performed using Gray 
Mapping in a quadrature phase shift keying device 220. 

30 Gray mapping is used to prevent that differential detection 
phase errors smaller than 135 degrees cause double bit er- 
rors at the receiver. 

Differential phase encoding of the K phases is performed in 
35 a differential phase encoder 222. At this stage of process- 
ing, the K phases phi [k] generated by the QPSK Gray mapper 
are differentially encoded. In principal, a feedback loop 
224 calculates a cumulative sum over all K phases. As start- 
ing point for the first computation (k = 0) the phase of the 
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reference carrier 226 is used. A switch 228 is provided in 
order to provide either the absolute phase of the reference 
subcarrier 226 or the phase information encoded onto the 
preceding (i.e. z 1 , where z 1 denotes the unit delay opera- 
5 tor) subcarrier to a summing point 230. At the output of the 
differential phase encoder 222, the phase information 
theta[k] with which the respective subcarriers are to be en- 
coded is provided. In preferred embodiments of the present 
invention, the subcarriers of a MCM symbol are equally 
10 spaced in the frequency axis direction. 

The output of the differential phase encoder 222 is con- 
nected to a unit 232 for generating complex subcarrier sym- 
bols using the phase information theta[k] . To this end, the 
15 K differentially encoded phases are converted to complex 
symbols by multiplication with 

factor * e j * [2 * pi * (theta[k]+PHl)] (Eq.l) 

20 wherein factor designates a scale factor and phi designates 
an additional angle. The scale factor and the additional an- 
gle PHI are optional. By choosing PHI = 45° a rotated DQPSK 
signal constellation can be obtained. 

25 Finally, assembly of a MCM symbol is effected in an assem- 
bling unit 234. One MCM symbol comprising N FFT subcarriers is 
assembled from N FFT — K— 1 guard band symbols which are "zero", 
one reference subcarrier symbol and K DQPSK subcarrier sym- 
bols. Thus, the assembled MCM symbol 200 is composed of K 

30 complex values containing the encoded information, two guard 
bands at both sides of the N FFT complex values and a refer- 
ence subcarrier symbol. 

The MCM symbol has been assembled in the frequency domain. 
35 For transformation into the time domain an inverse discrete 
Fourier transform (idft) of the output of the assembling 
unit 234 is performed by a transformator 236. In preferred 
embodiments of the present invention, the transformator 236 
is adapted to perform a fast Fourier transform (FFT) . 
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Further processing of the MCM signal in the transmitter as 
well as in the receiver is as described above referring to 
Figure 1. 

5 

At the receiver a de-mapping device 142 (Figure 1) is needed 
to reverse the operations of the mapping device described 
above referring to Figure 3. The implementation of the de- 
mapping device is straightforward and, therefore, need not 
10 be described herein in detail. 

The differential mapping along the frequency axis direction 
is suitable for multi -carrier (ofcm) digital broadcasting 
over rapidly changing multi path channels, in accordance 

15 with this mapping scheme, there is no need for a channel 
stationarity exceeding one multi -carrier symbol. However, 
differential mapping into frequency axis direction may cre- 
ate a new problem, in multi path environments, path echoes 
succeeding or preceding the main path can lead to systematic 

20 phase offsets between sub-carriers in the same MCM symbol. 
Thus, it will be preferred to provide a correction unit in 
order to eliminate such phase offsets. Because the channel 
induced phase offsets between differential demodulated sym- 
bols are systematic errors, they can be corrected by an al- 

25 gorithm. In principle, such an algorithm must calculate the 
echo induced phase offset from the signal space constella- 
tion following the differential demodulation and subse- 
quently correct this phase offset. 

30 Examples for such echo phase correction algorithms are de- 
scribed at the end of this specification referring to Fig- 
ures 12 to 15. 

Next, the fine frequency synchronization in accordance with 
35 the present invention will be described referring to Figures 
4 to 8. As mentioned above, the fine frequency synchroniza- 
tion in accordance with the present invention is performed 
after completion of the coarse frequency synchronization. 
Preferred embodiments of the coarse frequency synchroniza- 




- 16 - 



tion which can be performed by the symbol frame/carrier fre- 
quency synchronization unit 134 are described hereinafter 
referring to Figures 9 and 10 after having described the 
fine frequency synchronization in accordance with the pres- 
5 ent invention. 



nl 



With the fine frequency synchronization in accordance with 
the present invention frequency offsets which are smaller 
than half the sub-carrier distance can be corrected. Since 

10 the frequency offsets are low and equal for all sub-carriers 
the problem of fine frequency synchronization is reduced to 
sub-carrier level. Figure 4 is a schematical view of MCM 
symbols 200 in the time-frequency plane. Each MCM symbol 200 
consists of 432 sub-carrier symbols C x to C 432 . The MCM sym- 

15 bols are arranged along the time axis, the first MCM symbol 
200 shown in Figure 4 having associated therewith a time T lf 
the next MCM symbol having associated therewith a time T 2 
and so on. In accordance with a preferred embodiment of the 
present invention, the fine frequency synchronization is 

20 based on a phase rotation which is derived from the same 
sub-carrier of two MCM symbols which are adjacent in the 
time axis direction, for example CjT x and CjT 2 . 

in the following, the present invention is described refer- 
25 ring to QPSK mapping (QPSK = Quadrature Phase Shift Keying). 
However, it is obvious that the present invention can be ap- 
plied to any MPSK mapping, wherein M designates the number 
of phase states used for encoding, for example 2, 4, 8, 16 

30 

Figure 5 represents a complex coordinate system showing a 
QPSK constellation for each sub-carrier with frequency off- 
set. The four possible phase positions of a first MCM sym- 
bol, MCM-symbol-1 are shown at 300. Changing from the sub- 
35 carrier (sub-carrier n) of this MCM symbol to the same sub- 
carrier of the next MCM symbol, MCM-symbol -2 , the position 
in the QPSK constellation will be unchanged in case there is 
no frequency offset, if a frequency offset is present, which 
is smaller than half the distance between sub-carriers, as 
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mentioned above, this frequency offset causes a phase rota- 
tion of the QPSK constellation of MCM-symbol-2 compared with 
MCM-symbol-1. The new QPSK constellation, that is the four 
possible phase positions for the subject sub-carrier of MCM- 
symbol-2 are shown at 302 in Figure 5- This phase rotation 0 
can be derived from the following equation: 



C n designates the QPSK constellation of a sub-carrier n in a 
mcm symbol, n is an index running from 1 to the number of 
active sub-carriers in the MCM symbol- Information regarding 

the frequency offset is contained in the term e^ 2 ^"^ of 
equation 2. This frequency offset is identical for all sub- 
carriers- Therefore, the phase rotation 8 is identical for 
all sub-carriers as well. Thus, averaging overall sub- 
carrier of a MCM symbol can be performed. 

Figure 6 shows a block diagram of a MCM receiver in which 
the present invention is implemented. An analog/digital con- 
verter 310 is provided in order to perform an analog/digital 
conversion of a down-converted signal received at the re- 
ceiver front end 132 (Figure 1). The output of the ana- 
log/digital converter 310 is applied to a low path filter 
and decimator unit 312. The low path filter is an impulse 
forming filter which is identical to an impulse forming fil- 
ter in the MCM transmitter. In the decimator, the signal is 
sampled at the MCM symbol frequency. As described above re- 
ferring to Figure 1, guard intervals in the MCM signal are 
removed by a guard interval removal unit 132. Guard inter- 
vals are inserted between two MCM symbols in the MCM trans- 
mitter in order to avoid intersymbol interference caused by 
channel memory. 

The output of the guard interval removal unit 132 is applied 
to a MCM demodulator 314 which corresponds to the fast Fou- 
rier transformator 140 shown in Figure 1. Following the MCM 




(Eq.2) 
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demodulator 314 a differential decoding unit 316 and a de- 
mapping unit 318 are provided. In the differential decoding 
unit 316, phase information is recovered using differential 
decoding, in the demapping unit 318, demapping along the 
frequency axis direction is performed in order to recon- 
struct a binary signal from the complex signal input into 
the demapping unit 318. 

The output of the Mofl demodulator 314 is also applied to 
fine frequency error detector 320. The fine frequency error 
detector 320 produces an frequency error signal from the 
output of the MCM demodulator . in the depicted embodiment, 
the output of the finje frequency error detector 320 is ap- 

controlled oscillator 322 via a loop 
filter 324 is a low pass filter for 
interference portions of a higher 
Ijbwly varying error signal. The numer- 
ator 322 produces a carrier signal on 
I zered error signal. The carrier signal 
produced by the numerical controlled oscillator 322 is used 
for a frequency correction which is performed by making use 
of a complex multiplier 326. The inputs to the complex mul- 
tiplier 326 are the output of the low pass filter and deci- 
mator unit 312 and :he output of the numerical controlled 
oscillator 322. 



plied to a numerical 
filter 324. The loop 
f i 1 teri ng supe ri mposefd 
frequency from the si 
cal controlled oscil 
the basis of the fil 
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A description of a preferred embodiment of the fine fre- 
quency error detector 320 is given hereinafter referring to 
Figure 7. 

The fine frequency error detfector 320 comprises a differen- 
tial detector in the time axis 330. The output of the MCM 



demodulator 314, i.e. the 
Transform) is applied to the 
tector 330 which performs a 
time axis in order to derive 
set from the same sub-car rie 



T output (fft = Fast Fourier 
input of the differential de- 
differential detection in the 
information on a frequency off- 
r of two subsequently arriving 



MCM symbols, in the embodiment shown in Figure 7, the number 
of active sub-carriers is 432. Thus, the differential detec- 
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tor 330 performs a /correlation between the first and the 
433 rd sample- The hirst sample is associated with MCM- 
symbol-1 (Figure 5) J whereas the 433rd sample is associated 
with MCM-symbol -2 (JFigure 5). However, both these samples 
are associated with the same sub-carrier. 
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To this end, the input signal Y k is applied to a z^-block 
332 and thereafter to a unit 334 in order to form the com- 
plex conjugate of the output of the z^-block 332. A complex 
multiplier 336 is provided in order to multiply the output 
of the unit 334 by the input signal Y k . The output of the 
multiplier 336 is a signal Z k . 



m 



15 



The function of the differential 
pressed as follows: 



detector 330 can be ex- 



Z k = Y k+K ' Yk 



(Eq.3) 



hah 
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Y = [Y x ,Y 2 ...,Y k ...] 

Y = [c x /T x ,C 2 / T x ,...,C 432 / T X ,C X I T 2 ...] 



(Eq.4) 



CEq.5) 



25 



Y k designates the output of the MCM modulator 314, i.e. the 
input to the differential detector 330, at a time k. z k des- 
ignates the output of the differential detector 330. K des- 
ignates the number of active carriers. 



30 



35 



The output z k of the differential detector 330 contains a M- 
fold uncertainty corresponding to codeable phase shifts, in 
case of the QPSK this Mi-fold uncertainty is a 4-fold uncer- 
tainty, i.e. 0°, 90°, 180° and 270°. This phase shift uncer- 
tainty is eliminated from z k making use of a m-psk decision 
device 340. Such decision devices are known in the art and, 

5 described here in detail. The out- 

ce 340 (a k )* represents the complex 
e phase shift decided by the deci- 
tput of the decision device 340 is 



therefore, have not to b 
put of the decision devi 
conjugate of the codeab 
si on device 340. This oi 



correlated with the 
by performing a cor 
342. 
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>utput of the differential detector 330 
>lex multiplication using a multiplier 



10 
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The output the multiplier 342 represents the phase offset 
for the respective sub-aarriers . This phase offsets for the 
respective sub-carriers /are averaged over one MCM symbol in 
an averaging unit 344 irt accordance with a preferred embodi- 
ment of the present invention. The output of the averaging 
units 344 represent th< output of the fine frequency error 
detector 320. 

/ 

The mathematical description for this procedure is as fol- 
lows: 



" offset 



2ttKT m 



arg- 



K 

71 = 1 



(Eq.6) 



in accordance with preferred embodiments of the present in- 
vention, the frequency control loop has a backward struc- 
20 ture. in the embodiment shown in Figure 6, the feedback loop 
is connected between the output of the MCM demodulator 314 
and the input of the guard interval removal unit 132. 

in Figure 8, a block diagram of a MCM receiver comprising a 
25 coarse frequency correction unit 350 and a fine frequency 
correction unit as described above is shown. As shown in 
Figure 8, a common complex multiplier 326 can be used in or- 
der to perform the coarse frequency correction and the fine 
frequency correction. As shown in Figure 8, the multiplier 
30 326 can be provided preceding the low pass filter and deci- 
mator unit 312. Depending on the position of the multiplier 
326, a hold unit has to be provided in the fine frequency 
synchronization feedback loop. In an alternative embodiment, 
it is possible to use two separate multipliers for the 
35 coarse frequency correction and for the fine frequency cor- 
rection. In such a case, the multiplier for the coarse fre- 
quency correction will be arranged preceding the low path 
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filter and decimator unit, whereas the multiplier for the 
fine frequency correction will be arranged following the low 
path filter and decimator unit. 

5 Following, preferred embodiments for implementing a coarse 
frequency synchronization will be described referring to 
Figures 9 and 10. 



10 



15 



20 




30 



35 



As it is shown in Figure 9, the output of the receiver front 
end 132 is connected to an analog/digital converter 310. The 
down-converted MCM signal is sampled at the output of the 
analog/digital converter 310 and is applied to a 
frame/timing synchronization unit 360. in a preferred em- 
bodiment, a fast running automatic gain control (AGC) (not 
shown) is provided preceding the frame/timing synchroniza- 
tion unit in order to eliminate fast channel fluctuations. 
The fast AGC is used in addition to the normally slow AGC in 
the signal path, in the case of transmission over a mul- 
tipath channel with long channel impulse response and fre- 
quency selective fading. The fast AGC adjusts the average 
amplitude range of the signal to the known average amplitude 
of the reference symbol . 

As described abovd, the frame/timing synchronization unit 
uses the amplitude-lmodulated sequence in the received signal 
in order to extract the framing information from the MCM 
signal and furthen to remove the guard intervals therefrom. 
After the frame/timing synchronization unit 360 it follows a 
coarse frequency synchronization unit 362 which estimates a 
coarse frequency dffset based on the amplitude-modulated se- 
quence of the reference symbol of the MCM signal, in the 



coarse frequency 
set of the carrie 
frequency in the 
form a frequency 
quency offset cor 
pi ex mul ti pi i cati 



synchronization unit 362, a frequency off- 
frequency with respect to the oscillator 
MCM receiver is determined in oder to per- 
offset correction in a block 364. This fre- 
rection in block 364 is performed by a corn- 
on . 
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The output of the frequency offset correction block 364 is 
applied to the MCM demodulator 366 formed by the Fast Fou- 
rier Transformator 140 and the carrier-bit mapper 142 shown 
in Figure 1. 

5 

In order to perform the coarse frequency synchronization de- 
scribed herein, an amplitude-demodulation has to be per- 
formed on a preprocessed MCM signal. The preprocessing may 
be, for example, the down-conversion and the analog/digital 
10 conversion of the MCM signal. The result of the amplitude- 
demodulation of the preprocessed MCM signal is an envelope 
representing the amplitude of the MCM signal. 

For the amplitude demodulation a simple alpha^ beta min . 
15 method can be used. This method is described for example in 
Palacherla A.: DSP-yP Routine Computes Magnitude, EDN, Octo- 
ber 26, 1989; and Adams, w. T. , and Bradley, J.: Magnitude 
Approximations for Microprocessor Implementation, IEEE Mi- 
cro, Vol. 3, No. 5, October 1983. 

20 

It is clear that amplitude determining methods different 
from the described alpha max+ beta min . method can be used. For 
simplification, it is possible to reduce the amplitude cal- 
culation to a detection as to whether the current amplitude 
25 is above or below the average amplitude. The output signal 
then consists of a -1/+1 sequence which can be used to de- 
termine a coarse frequency offset by performing a correla- 
tion. This correlation can easily be performed using a sim- 
ple integrated circuit (IC) . 

30 

in addition, an oversampling of the signal received at the 
RF front end can be performed. For example, the received 
signal can be expressed with two times oversampling. 

35 in accordance with a first embodiment, a carrier frequency 
offset of the MCM signal from an oscillator frequency in the 
MCM receiver is determined by correlating the envelope ob- 
tained by performing the amplitude-demodulation as described 
above with a predetermined reference pattern. 
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10 



in case there is no frequency offset, the received reference 
symbol r(k) will be: 



r(k) = S^Ck) + n(k) 



(Eq.7) 



wherein n(k) designates "additive Gaussian noise" and de- 
notes the AM sequence which has been sent, in order to sim- 
plify the calculation the additive Gaussian noise can be ne- 
glected- It follows: 



r(k)^S AM (k) 



(Eq.8) 



In case a constant frequency offset Af is present, the re- 
15 ceived signal will be: 



CEq.9) 



information regarding the frequency offset is derived from 
20 the correlation of the received signal r(k) with the AM se- 
quence which is known in the receiver: 



25 



Ar=l *=1 

Thus, the frequency offset is: 



(Eq.10) 



f L 
2 



2aT, 



arg 



MCM 



*=1 



f L 



2nT h 



-arg 



MCM 



jt=i 



(Eq.ll) 



30 Since the argument of IS^Ck)! 2 is zero the frequency offset 
is: 
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1 



2 



arg £F(A:).S; 



(Eq.12) 



k=\ 

V 



in accordance with a second embodiment of the coarse fre- 
quency synchronization algorithm, a reference symbol com- 
prising at least two identical sequences 370 as shown in 
Figure 10 is used. Figure 10 shows the reference symbol of a 
MCM signal having two identical sequences 370 of a length of 
L/2 each. L designates the number of values of the two se- 
quences 370 of the reference symbol . 

As shown in Figure 10, within the amplitude-modulated se- 
quence, there are at least two identical sections devoted to 
the coarse frequency synchronization. Two such sections, 
each containing L/2 samples, are shown at the end of the am- 
plitude-modulated sequence in Figure 10. The amplitude- 
modulated sequence contains a large number of samples. For a 
non-ambiguous observation of the phase, only enough samples 
to contain a phase rotation of 2n should be used. This num- 
ber is defined as L/2 in Figure 10. 

Following, a mathematical derivation of the determination of 
a carrier frequency deviation is presented, in accordance 
with Figure 10, the following equation applies for the two 
identical sequences 370: 



If no frequency offset is present, the following equation 14 
will be met by the received signal: 




(Eq.13) 





(Eq.14) 
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r(k) designates the values of the identical sequences 
an index from one to L/2 for the respective samples. 



k is 



if there is a frequency offset of, for example, Af, the re- 
5 ceived signal is: 



r(k) = r(k)e J2 ** kTua ' 



CEq.15) 



10 



15 



r(k + -) = r(k)e 



(Eq.16) 



rCk) designates sample values of the received portion which 
are based on the identical sequences. Information regarding 
the frequency offset is derived from the correlation of the 
received signal f (k + l / 2) with the received signal r(k). 
This correlation is given by the following equation: 



k=\ V Zy k=\ 



CEq-17) 



f* designates the complex conjugate of the sample values of 
20 the portion mentioned above. 

Thus, the frequency offset is 



25 



f L 



Af = 



2k - T, 



arg 



MCM 



f L 
2 



2t ^ T MCM 



arg 



L 

(Eq.18) 



Since the argument of \r(kf equals zero, the frequency offset 
becomes 
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f L 



Af = 



2ji - r M 



arg 



r\k + 



(Eq.19) 



Thus, it is clear that in both embodiments, described above, 
the frequency position of the maximum of the resulting out- 
put of the correlation determines the estimated value of the 
offset carrier. Furthermore, as it is also shown in Figure 
9, the correction is performed in a feed forward structure. 



case of a channdl with strong reflections, for example 



10 1 due to a high build 



third embodiment of 
values of the two 

15 dance with a second 
sponding values of 
corresponding to sai 
ence symbol. This wes 
correctly determinin 

20 description of this 



ling density, the correlations described 
above might be insufficient for obtaining a suitable coarse 
frequency synchronisation. Therefore, in accordance with a 

the present invention, corresponding 
portions which are correlated in accor- 
embodiment, can be weighting with cor re- 
stored predetermined reference patterns 
d two identical sequences of the refer- 
iighting can maximize the probability of 
the frequency offset. The mathematical 
i/eighting is as follows: 



25 



Af = 



2k - 



arg 



( L 

2 

I 

Jt = l _ 



(Eq.20) 

designates the amplitude-modulated sequence which is 
known in the receiver, and s* m designates the complex conju- 
gate thereof. 



* AM 



if the above correlations are calculated in the frequency 
30 domain, the amount of 



(Eq.21) 
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is used rather than the argument. This amount is maximized 
as a function of a frequency correction. The position of the 
maximum determines the estimation of the frequency devia- 
5 tion. As mentioned above, the correction is performed in a 
feed forward structure. 



10 



Preferred embodiments for performing an echo phase offset 
correction when using a differential mapping in the fre- 
quency axis will be described hereinafter referring to Fig- 
ures 12 to 15. 




Systematic phase 
envi ronments ma> 
symbol . Thi s ph 
lating the MCM s 
to make use of 
shifts stemming 



shifts stemming from echoes in multipath 
occur between subcarriers in the same MCM 
offsets can cause bit errors when demodu- 
ymbol at the receiver. Thus, it is preferred 
algorithm to correct the systematic phase 
from echoes in multipath environments. 




ase 



In Figure 12, scatter 



diagrams at the output of a differen- 



25 



tial demapper of a MCM receiver are shown. As can be seen 
from the left part of Figure 12, systematic phase shifts be- 
tween subcarriers in :he same MCM symbol cause a rotation of 
the demodulated phase Ishifts with respect to the axis of the 
complex coordinate system. In the right part of Figure 12, 
the demodulated phase Ishifts after having performed an echo 
phase offset correction are depicted. Now, the positions of 



30 



the signal points are 
pi ex coordinate system 
modulated phase shifts 
ti vely . 



substantially on the axis of the com- 
These positions correspond to the 
of 0°, 90°, 180° and 270°, respec- 



35 




An echo phase offset correction algorithm (EPOC algorithm) 
must calculate the echo induced phase offset from the signal 
space constellation following the differential demodulation 
and subsequently correct this phase offset. 

For illustration purposes, one may think of the simplest al- 
gorithm possible which eliminates the symbol phase before 



* 
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computing the mean of all phases of the subcarriers. To il- 
lustrate the effect/ of such an EPOC algorithm, reference is 
made to the two sca/tter diagrams of subcarriers symbols con- 
tained in one MCM symbol in Figure 12. This scatter diagrams 
5 have been obtained/ as result of an MCM simulation. For the 
simulation a channel has been used which might typically 
show up in single frequency networks. The echoes of this 
channel stretched/to the limits of the MCM guard interval. 
The guard interval was chosen to be 2 5% of the MCM symbol 
10 duration in this case. 



Figure 13 represents a block diagram for illustrating the 
( position and the functional iry of an echo phase offset cor- 
15 rection device in a MCM receiver. The signal of a MCM trans- 
mitter is transmitted through the channel 122 (Figures 1 and 
13) and received at the receiver frontend 132 of the MCM re- 
ceiver. The signal processi/ng between the receiver frontend 
and the fast Fourier tranOTormator 140 has been omitted in 
20 Figure 13. The output of ithe fast Fourier transformator is 
applied to the de-mapper, /which performs a differential de- 
mapping along the frequency axis. The output of the de- 
mapper are the respective phase shifts for the subcarriers. 
The phase offsets of thjs phase shifts which are caused by 
25 echoes in multipath environments are visualized by a block 
400 in Figure 13 which shows an example of a scatter diagram 
of the subcarrier symbojls without an echo phase offset cor- 
rection. 

30 The output of the de-mapper 142 is applied to the input of 
an echo phase offset correction device 402. The echo phase 
offset correction device 402 uses an EPOC algorithm in order 
to eliminate echo phase offsets in the output of the de- 
mapper 142. The result is shown in block 404 of Figure 13, 

35 i.e. only the encoded phase shifts, 0°, 90°, 180° or 270° 
are present at the output of the correction device 402. The 
output of the correction device 402 forms the signal for the 
metric calculation which is performed in order to recover 
the bitstream representing the transmitted information. 



# # 
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}/ w /a first embodiment of an EPOC algorithm and a device for 
I performing same ps now described referring to Figure 14. 

5 The first embodiment of an EPOC algorithm starts from the 
assumption that every received differentially decoded com- 
plex symbol is rotated by an angle due to echoes in the mul- 
tipath channel. For the subcarriers equal spacing in fre- 
quency is assumed since this represents a preferred embodi- 
10 ment. If the subcarriers were not equally spaced in fre- 
quency, a correction factor would have to be introduced into 
the EPOC algorithm. 

Figure 14 shows the correction device 402 (Figure 13) for 
15 performing the first embodiment of an EPOC algorithm. 



From the output of the de-mapper 142 which contains an echo 
phase offset as shown for example in the left part of Figure 
12, the phase shifts related to transmitted information must 

20 first be discarded. To this end, the output of the de-mapper 
142 is applied to a discarding unit 500. in case of a DQPSK 
mapping, the discarding unit can perform a " (.) 4M operation. 
The unit 500 projects all received symbols into the first 
quadrant. Therefore, the phase shifts related to transmitted 

2 5 information is eliminated from the phase shifts representing 
the subcarrier symbols. The same effect could be reached 
with a modulo-4 operation. 

Having eliminated the information related symbol phases in 
30 unit 500, the first approach to obtain an estimation would 
be to simply compute the mean value over all symbol phases 
of one MCM symbol. However, it is preferred to perform a 
threshold decision before determining the mean value over 
all symbol phases of one MCM symbol. Due to Rayleigh fading 
35 some of the received symbols may contribute unreliable in- 
formation to the determination of the echo phase offset. 
Therefore, depending on the absolute value of a symbol, a 
threshold decision is performed in order to determine 
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whether the symbol should contribute to the estimate of the 
phase offset or not. 

Thus, in the embodiment shown in Figure 14, a threshold de- 
5 cision unit 510 is included. Following the unit 500 the ab- 
solute value and the argument of a differentially decoded 
symbol is computed in respective computing units 512 and 
514. Depending on the absolute value of a respective symbol, 
a control signal is derived. This control signal is compared 

10 with a threshold value in a decision circuit 516. if the ab- 
solute value, i.e. the control signal thereof, is smaller 
than a certain threshold, the decision circuit 516 replaces 
the angle value going into the averaging operation by a 
value equal to zero. To this end, a switch is provided in 

15 order to disconnect the output of the argument computing 
unit 514 from the input of the further processing stage and 
connects the input of the further processing stage with a 
unit 518 providing a constant output of "zero". 

20 An averaging unit 520 is provided in order to calculate a 
mean value based on the phase offsets <p. determined for the 
individual subcarrier symbols of a MCM symbol as follows: 

W = i/*I>i (Eq.22) 

i=l 

25 

in the averaging unit 520, summation over K summands is per- 
formed. The output of the averaging unit 520 is provided to 
a hold unit 522 which holds the output of the averaging unit 
520 K times. The output of the hold unit 522 is connected 
30 with a phase rotation unit 524 which performs the correction 
of the phase offsets of the K complex signal points on the 
basis of the mean value ^. 

The phase rotation unit 524 performs the correction of the 
35 phase offsets by making use of the following equation: 



v' k = v k • e 



(Eq.23) 
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in this equation, v' k designates the K phase corrected dif- 
ferentially decoded symbols for input into the soft-metric 
calculation, whereas v k designates the input symbols. As 
5 long as a channel which is quasi stationary during the dura- 
tion of one MCM symbols can be assumed, using the mean value 
over all subcarriers of one MCM symbol will provide correct 
results. 

10 A buffer unit 527 may be provided in order to buffer the 
complex signal points until the mean value of the phase off- 
sets for one MCM symbol is determined. The output of the 
phase rotation unit 524 is applied to the further processing 
stage 526 for performing the soft-metric calculation. 

15 

With respect to the results of the above echo phase offset 
correction, reference is made again to Figure 12. The two 
plots stem from a simulation which included the first em- 
bodiment of an echo phase offset correction algorithm de- 

20 scribed above. At the instant of the scatter diagram snap- 
shot shown in the left part of Figure 12, the channel obvi- 
ously distorted the constellation in such a way, that a sim- 
ple angle rotation is a valid assumption. As shown in the 
right part of Figure 12, the signal constellation can be ro- 

25 tated back to the axis by applying the determined mean value 
for the rotation of the differentially detected symbols. 

A second embodiment of an echo phase offset correction algo- 
rithm is described hereinafter. This second embodiment can 
30 be preferably used in connection with multipath channels 
that have up to two strong path echoes. The algorithm of the 
second embodiment is more complex than the algorithm of the 
first embodiment. 

35 what follows is a mathematical derivation of the second em- 
bodiment of a method for echo phase offset correction. The 
following assumptions can be made in order to ease the ex- 
planation of the second embodiment of an EPOC algorithm. 
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In this embodiment, the guard interval of the MCM signal is 
assumed to be at least as long as the impulse response h[q], 
q = 0, 1, . .., Qh-1 of the multipath channel. 



At the transmitter every MCM symbol is assembled using fre- 
quency axis mapping explained above. The symbol of the ref- 
erence subcarrier equals 1, i.e. 0 degree phase shift. The 
optional phase shift PHI equals zero, i.e. the DQPSK signal 
10 constellation is not rotated. 

Using an equation this can be expressed as 



a k - a k-i a k 



inc 



15 (Eq.24) 



with 
k 



index k = 1,2 
rier ; 



K of the active subcar- 



inc 




complex phase increment symbol; m=0,l,2,3 
is the QPSK symbol number which is derived 
from Gray encoding pairs of 2 Bits; 



aO = 1 



symbol of the reference subcarrier. 



25 



At the DFT output of the receiver the decision variables 



(Eq.25) 



30 are obtained with 




CEq.26) 



i = 0 



35 



being the DFT of the channel impulse response h[q] at posi- 
tion k. 



# • 
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With | ak | 2 = 1 the differential demodulation yields 

n = e k • e;_ x = ai nc H k Hl^ (Eq.27) 

5 For the receiver an additional phase term <p k is introduced, 
which shall be used to correct the systematic phase offset 
caused by the channel. Therefore, the final decision vari- 
able at the receiver is 

10 v k = v k . e J * = a? c • - H k • 

(Eq.28) 

As can be seen from the Equation 28, the useful information 
a k nc is weighted with the product e j<Pk • H k • (rotation and 

15 effective transfer function of the channel). This product 
must be real-valued for an error free detection. Considering 
this, it is best to choose the rotation angle to equal the 

negative argument of H k - H* k _ im To derive the desired algo- 
rithm for 2-path channels, the nature of H k * is investi- 
20 gated in the next section. 

It is assumed that the 2-path channel exhibits two echoes 
with energy content unequal zero, i.e. at least two dominant 
echoes. This assumption yields the impulse response 

25 

%] = cMq] + c 2 S 0 [q - g 0 ] (Eq.29) 
with 

c lf c 2 : complex coefficients representing the path 

30 echoes; 

q 0 : delay of the second path echo with respect 

to the first path echo; 
£ 0 : Dirac pulse; S Q [k] = 1 for k = 0 

S 0 [k] =0 else 



35 
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The channel transfer function is obtained by applying a DFT 
to Equation 29: 



H k = H\e * 



= c, + c 0 



-J — kq 0 



(Eq.30) 



With Equation 30 the effective transfer function for differ- 
ential demodulation along the frequency axis is: 



10 



c ± + c 2 e K 
V J 



c x + c 2 e 



= C a + C K CO 



s[| q 0 (2 k - 1) 



(Eq.31) 



15 



Assuming a noise free 2-path channel, it can be observed 
from Equation 31 that the symbols on the receiver side are 
located on a straight line in case the symbol 1+jO has been 
send (see above assumption). This straight line can be char- 
acterized by a point 



2/r 
-3 — Go 



c a = N + |c 2 | • e * 



(Eq.32) 



20 and the vector 



c b 20-^2 



-J — <7o 



(Eq.33) 



which determines its direction. 

25 

with the above assumptions, the following geometric deriva- 
tion can be performed. A more suitable notation for the geo- 
metric derivation of the second embodiment of an EPOC algo- 
rithm is obtained if the real part of the complex plane is 
30 designated as x = Re{z}, the imaginary part as y = lm{z} , 
respectively, i.e. z = x+jy. with this new notation, the 
straight line, on which the received symbols will lie in 
case of a noise-free two-path channel, is 
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f (x) = a + b - x 



CEq.34) 



wi th 



a = Im{cJ - 



Re{c a } 
Re{cJ 



(Eq.35) 



and 



10 



b = - 



Re{c a } - 




(Eq.36) 



Additional noise will spread the symbols around the straight 
line given by Equations 34 to 36. in this case Equation 36 
is the regression curve for the cluster of symbols. 



For the geometric derivation of the second embodiment of an 
EPOC algorithm, the angle <p k from Equation 28 is chosen to 

be a function of the square distance of the considered sym- 
bol from the origin: 



Equation 37 shows that the complete signal space is dis- 
torted (torsion), however, with the distances from the ori- 
25 gin being preserved. 

For the derivation of the algorithm of the second embodi- 
ment, f K {) has to be determined such that all decision vari- 
ables v' k (assuming no noise) will come to lie on the real 
30 axi s : 



15 



20 




(Eq.37) 



lm|[x + jf(x)) ■ e jf "^} = 0 



(Eq.38) 
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Further transformations of Equation 38 lead to a quadratic 
equation which has to be solved to obtain the solution for 

<P* ■ 

5 in case of a two-path channel , the echo phase offset correc- 
tion for a given decision variable v k is 



(Eq.39) 



10 wi th 



<P k = 



— a tan 



a + byl\v k \ 2 (l + b 2 ) - a 2 



- ab + j\v k \ z (l + b 2 ) - a 2 



for \v k \ 2 > 



a tan 



for 



\ v k\ < 



i + & 

a 2 
1 + b 2 

(Eq.40) 



15 From the two possible solutions of the quadratic equation 
mentioned above, Equation 40 is the one solution that cannot 
cause an additional phase shift of 180 degrees. 

) (The two plots in Figure 15 show the projection of the EPOC 

20 algorithm of the second embodiment for one quadrant of the 
complex plane. Depncted here is the quadratic grid in the 
sector |arg(z)| <W/4 and the straight line 
y = f(x) = a + jb • x irith a = -1.0 and b = 0.5 (dotted line), 
in case of a noise-rFree channel, all received symbols will 

25 lie on this straight! line if l+j0 was send. The circle shown 
in the plots determiVies the boarder line for the two cases 
of Equation 40. in tlte left part, Figure 15 shows the situa- 
tion before the projection, in the right part, Figure 15 
shows the situation after applying the projection algorithm. 

30 By looking on the left part, one can see, that the straight 
line now lies on the real axis with 2+j0 being the fix point 
of the projection. Therefore, it can be concluded that the 
echo phase offset correction algorithm according to the sec- 
ond embodiment fulfillsythe design goal. 
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Before the second embodiment of an EPOC algorithm can be ap- 
plied, the approximation line through the received symbols 
has to be determined, i.e. the parameters a and b must be 
5 estimated- For this purpose, it is assumed that the received 
symbols lie in sector |arg(z)| < n / 4 , if 1+jO was sent. If 

symbols other than 1+jO have been sent, a modulo operation 
can be applied to project all symbols into the desired sec- 
tor. Proceeding like this prevents the necessity of deciding 
10 on the symbols in an early stage and enables averaging over 
all signal points of one MCM symbol (instead of averaging 
over only H of all signal points). 

For the following computation rule for the EPOC algorithm of 
15 the second embodiment, x i is used to denote the real part of 
the i-th signal point and y., for its imaginary part, respec- 
tively (i = 1, 2,..., K) . Altogether, K values are available 
for the determination. By choosing the method of least 
squares, the straight line which has to be determined can be 
20 obtained by minimizing 



The solution for Equation 41 can be found in the laid open 
25 literature. It is 



(a, jb) = arq min jT (y. - (a + B • xjf 



(Eq.41) 



K 



b = 



i=l 



a 



= y — x - b 



(Eq.42) 



K 



1=1 



with mean values 



30 




(Eq.43) 
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If necessary, an estimation method with higher robustness 
can be applied. However, the trade-off will be a much higher 
computati onal compl exi ty . 

5 To avoid problems with the range in which the projection is 
applicable, the determination of the straight line should be 
separated into two parts. First, the cluster's centers of 
gravity are moved onto the axes, following, the signal space 
is distorted. Assuming that a and b are the original parame- 
10 ters of the straight line and a is the rotation angle, f K () 
has to be applied with the transformed parameters 



Besides the two EPOC algorithms explained above section, 
different algorithms can be designed that will, however, 
most likely exhibit a higher degree of computational com- 
plexity. 




(Eq.44) 
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ABSTRACT 



method and an apparatus relate to a fine frequency syn- 
chronization compensating for a carrier frequency deviation 
from an oscillator requency in a multi-carrier demodulation 
system of the type capable of carrying out a differential 
phase decoding of mjlti -carrier modulated signals, the sig- 
nals comprising a plurality of symbols, each symbol being 
defined by phase differences between simultaneous carriers 
having different frequencies. A phase difference between 



phases 
mined. 



of the same 
Thereafter, 



carrier in different symbols is deter- 
a frequency offset is determined by 
eliminating phase shift uncertainties related to the trans- 
mitted information from the phase difference making use of a 
M-PSK decision device. Finally, a feedback correction of the 
carrier frequency deviation is performed based on the deter- 
mined frequency off sat. Alternatively, an averaged frequency 
offset can be determined by averaging determined frequency 
offsets of a plurality of carriers. Then, the feedback cor- 
rection of the frequency deviation is performed based on the 
averaged frequency offset. 



